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Abstract 

We study the cooling of a dielectric nanoscale particle trapped in an optical cavity. We derive the frictional force for motion 
in the cavity field, and show that the cooling rate is proportional to the square of oscillation amplitude and frequency. Both the 
radial and axial centre-of-mass motion of the trapped particle, which are coupled by the cavity field, are cooled. This motion 
is analogous to two coupled but damped pendulums. Our simulations show that the nanosphere can be cooled to e _1 of its 
initial momentum over time scales of hundredths of milliseconds. 
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The study of micro- and nanomechanical oscillators [U EJ [3j [4] , and particularly their quantum mechanical motion 
[SJEIIZ], is a rapidly developing field which promises insights into the boundary between quantum and classical worlds. 
Also, their sensitivity to the environment appears promising for the development of chemical sensors with single atom 
sensitivity [8]. Cavity opto- mechanics is an important area within this field, in which the mechanical motion of at 
least one degree of the freedom of an oscillator is damped or cooled by interaction with the field of an optical cavity 
[H EJ [lOj [TTJ [12] . Cooling of the mechanical motion works by coupling an optical cavity to the oscillator so that it 
selectively scatters blue shifted photons of a probe beam out of the cavity with respect to the incident photons. By 
conservation of energy, the mechanical energy of the oscillator-cavity system is reduced. A range of optomechanical 
oscillators and cooling schemes have been now been realized experimentally [9 , and currently there is a strong impetus 
towards reaching the quantum regime, where only a few of the quantised states of at least one degree of freedom are 
occupied. The cooling of atomic species using laser cooling is now well established with the creation of nanoKelvin 
temperature gases, which has led to the realization of quantum degeneracy in gases. For molecular and atomic species 
that cannot be laser cooled, cavity cooling of a trapped species appears attractive because it does not rely on the 
detailed internal level structure p~3j [14] . This scheme has already been used to cool a trapped atom [15] , an ion [16] , 
and atomic gases [T7] . 




FIG. 1: Schematic diagram for cooling a small dielectric nanoscale particle in a high finesse cavity. 



I. INTRODUCTION 



In this paper we consider the cavity cooling of a nanoscale particle trapped at the center of a high finesse Fabry- 
Perot cavity in a vacuum. Here a center-of-mass oscillation occurs due to trapping in the periodic potential of the 
interference pattern created inside the cavity, while cooling occurs by interaction with the field mediated by the cavity. 
Unlike many opto-mechanical cooling schemes that utilize radiation pressure, this scheme uses the dipole force. Cavity 
cooling in this way is attractive because a trapped particle can be effectively isolated from the environment. This is 
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unlike the cooling of many cavity optomechanical schemes such as cantilever [4] or membranes [18] , which are directly 
physically connected to a large heat bath. 

We consider cavity cooling of a large polarizable particle that is trapped in the intra-cavity field of a high finesse 
cavity by modifying the model developed for a single atom in a 1-D cavity [19]. Because we only consider a single 
trapped particle, we are able to derive a velocity dependent frictional force for small oscillations around the antinode 
of the intra-cavity field. We then consider cooling in a realistic 2-D cavity which includes the damping of the axial 
and radial motion. We show that both degrees of freedom are coupled by the cavity field which acts to damp them. 
In addition to dipole force due to the cavity field, which acts to trap the particle, the nanoparticle is also subject to 
gravity, which is chosen to be along the cavity axis as shown in figure 1. 



II. 1-D COOLING OF A NANOSPHERE 



To explore cavity cooling we assume a field of amplitude ^ ext and frequency uo p incident on one of the two high 
finesse cavity mirrors of equal reflectivity, R, with conductivity a. The line width of the cavity is n and the 1/e 
lifetime of a photon in the cavity is ^ . We assume that the light is impedance matched to the cavity, that the cavity 
is stabilized such that only the TEMqo mode propagates, and that the nanoparticle is trapped at its beam waist. We 
first determine the field within the cavity in the presence of the single dielectric particle from the one-dimensional 
wave equation given by 

d 2 E(x,t) dE(x,t) 2 d 2 E(x,t) _ 

i d 2 P{x,t) ext dE(x,ty** 

7 dt 2 + 2/ " dt (1) 

where k = J- and n ext = ^ ext . 

Here E(x,t) and P(x,t) represent the sum of all possible allowed electric field and polarization modes for the 
cavity respectively. We assume that the cavity only operates on one of these modes where the electric field is 
E m (x, t) = ^(t)(e~ luJpt + c.c)cos(kx) and the polarization is given by P m = p(t)e~ luJpt + c.c. in the cavity. The external 
field which couples to this cavity mode is given by E™ 1 — £ ext (t}e~ lu}pt + c.c. We find the equation for the evolution 
of field amplitude of this mode by utilizing the orthogonality of cavity modes and multiplying equation (1) by E m (x, t) 
and averaging over the cavity volume V . 

Under the slowly varying envelope approximation, we obtain the 1-D equation of motion for the amplitude of the 
single mode field as 



-± = -[K-i(U (x) + A)]£ + K ext iext (2) 

where U (x) = aUp ™y kx is the position-dependent shift in cavity frequency induced by the polarizable particle and 
A = uo p — uo c is the cavity detuning from resonance. 

The equation of motion for the momentum is simply the dipole force that results from the intra-cavity field acting 
against gravity 

dP 

—j^- = — a\^\ 2 ksin2kx — mg (3) 

and the particle position is given by 



Tt = P * /m (4) 

It is not clear from the coupled equations ([3| and Q that there is a dissipative force on the trapped particle 
in the intracavity field. To understand the damping or heating of the trapped nanoparticle we derive an equation of 
motion with a velocity dependent frictional force. As we are only interested in trapped motion, we need only consider 
small amplitude oscillations when kx << 1. Additionally, because the center-of-mass motion is much smaller than 
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the round trip time for light in the cavity, we average equations ([2]), (|3| and Q over a time At from t — At to t. A 
suitable At is the cavity decay time 2tt / ' n. The averaged field £ over 2i\ j k is given by 

„ext text 

— (5) 



where the position dependent expression for the frequency shift U(x) is expanded around x = so that U(x) « 
^ p (i-(fex) 2 ) and ^ _ A ^ = ^ _ 2vir / KX (t), where i; = §. 

This enables us to find a time-averaged equation of motion for momentum. This is given by ^ = M^f « 

—2ak 2 ^ x, which leads to an equation of motion of the trapped nanosphere which is a type of Lienard differential 
equation [20]. This is given by 

= -n 2 (x + px 3 - ^x 2 v) - mg (6) 
at 1 k 

where Q — \/ 2a ^ 2 , S = J* f eJ„, 9 and /? = T ^ r P 9 , 7 A + SX! 91 . This equation differs from conventional oscillators 

such as the van der Pol or Duffing type in that the friction term, -x 2 , is position dependent. When /3 < the 

oscillatory motion of the sphere is damped and therefore the cavity detuning is A < — aco p /eoV. Optimal damping 
or cooling of the sphere occurs when 

A op -a(jp/e V - k. (7) 

An approximate expression for the energy loss rate averaged over an oscillation period 2tt/Q can also be determined. 
The rate of change in the total energy, E, of the particle is ^ = Mv^+M£} 2 xj^+MfiQ 2 x 3 , while the power can be 

determined by multiplying the equation of motion ^ by Mv to give Mv^- = —MQ 2 xv — MQ 2 Px 3 v+ 2 ^Mfl 2 v 2 x 2 . 
Equating equal terms gives the rate of change in total energy of the system as ^ = 2 -^-Mx 2 Q?v 2 . When averaged 
over an oscillation period, and assuming the damping is slow enough not to change the amplitude of oscillations during 
a cycle, we determine a time averaged exponential damping rate of 

r = ^xl (8) 

where xo is the amplitude of the oscillation. Because the amplitude will decrease as the particle cools at a fixed 
input intensity, so the damping rate also decreases, as has previously been noted [2T] , However, the cooling rate is 
proportional to the square of the oscillation frequency. 

To study the cooling of nanoscale polarizable particles, we consider a Si02 nanosphere of radius r = 100 nm and 
refractive index of n=1.45 trapped in an antinode of the intra-cavity field. A sphere of this size within the Rayleigh 
regime when r « A and can be treated as a dipole scatterer with polarizability a = 47reo ^ 2 ~j, r 3 , determined to be 
2.98 x 10 -32 Cm 2 V -1 . Its mass, determined from its density (p = 2198 Kg m -3 ), is 9.2 x 10 -18 kg. We consider a 
cavity of length L = 10 cm, with mirrors of reflectivity R = 0.99995 resulting in k = 4.71 x 10 5 rad/s. The intracavity 
field has a wavelength of 1064 nm with an input intensity of 100 mW and a beam radius of 50 microns in the cavity. 
The axial oscillation frequency for these parameters is 1.55 x 10 6 rad/s. The optimal detuning of the cavity from 
resonance is calculated from the above parameters to be —5.13 x 10 5 rad/s for cooling and 4.29 x 10 5 rad/s for heating. 
It is instructive to plot the energy damping/heating rate as a function of both detuning and position, as shown in 
figure 2. Peak negative values correspond to maximal cooling and peak positive values for maximal heating. Figure 
2 also shows that the damping/heating is larger when the particle is further from the anti-node of the interference 
pattern of the intracavity field, which has important consequences for cooling because it implies that the best damping 
will occur when the particle is not strongly confined and will undergo larger oscillations. 

To illustrate the damped oscillatory motion of the trapped nanosphere, we plot its velocity and position for an 
initial velocity of v = 3 cm s _1 , corresponding to an initial kinetic energy of E/k^ = 300if and x = 100 nm from 
trap center with optimal cavity detuning for damping. The equations of motion were integrated in time using a 
fourth-order Runge-Kutta method. The damped oscillatory motion is clearly shown, even for this short time period 
t = 1ms, and there is good agreement between the simplified time-averaged equation of motion ([6| and the exact 
model of equations ( 2|3|4 ). The oscillation frequency, largely determined by the input intensity, has been chosen to 



be smaller than the cavity line width cavity. However, an increase in intensity and thus oscillation frequency can be 



4 




I rtHfi 



FIG. 2: A 3-D plot showing the variation of energy damping rate as a function of cavity detuning A and oscillation amplitude 
of nanosphere position xq. Maximal cooling occurs at negative detuning. 
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FIG. 3: A phase space plot showing damped motion in position over 1 ms. Note that the cooling rate decreases as the particle 
becomes localized in the lattice potential. 



used to cool the particle in the resolved-sideband limit, of importance for cooling to the quantum ground state [22] . 
Figure 3 shows a phase space plot of the velocity and position trajectory of the particle trapped in a single anti-node. 
The inward spiral indicates damping with a much tighter spiral as the particle approaches x = 0, indicating a weaker 
frictional force as the particle approaches the ant i- node. It has been shown previously through simulations that atoms 
can be cooled more rapidly in cavity cooling schemes if the laser input intensity is lowered as the particle is cooled 
[2T] . This can now be seen to be due to the position dependence of the damping rate. For high intensities the particle 
is tightly confined and the amplitude of the oscillations are small resulting in a lower cooling rate. To illustrate the 
differences in cooling time for a constant intensity field, we plot the amplitude of the oscillating kinetic energy of 
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the particle as a function of time for a constant intensity of 1 mW input into the cavity and for an exponentially 
decaying input field, initially at 1 mW, with a time constant of 50 ms. We note that the lowest intensity that can be 
used is limited by gravity since eventually the well depth of the intra-cavity field becomes less than the gravitational 
potential. This occurs in figure 4 at approximately 0.35 s. Here, the cavity intensity must be held constant and the 
cooling rate where the cooling rate will be slower. 
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FIG. 4: Plot of the amplitude of the oscillatory energy as a function of time to compare of cooling times for a constant input 
intensity into the cavity of 1 mW (blue line) and an exponentially decreasing input intensity with e _1 decay constant of 50 ms 
(red line). 



III. COOLING IN 2-D 

We now consider the experimentally more realistic situation of cooling in a cavity which has field variations in the 
axial and in the radial direction. We now must also consider motion in the radial direction as well. For this case we 
assume the field input into the cavity is that of the lowest order transverse electromagnetic mode 



E(r,x,t) = A(t) 



(9) 



where ljq 



T. V 



r] = arctan — and uj(x) — c^owl + ( — ) 2 . We require that the nanosphere is 



trapped at the center of the cavity where the intra-cavity intensity is highest. At this position (xq = L/2) uj(x) « y/2i 



R(x) ~ 2xq and r] 



It is useful to recast this in trigonometric form as E(r,x,t) = | exp(— )[cos /ex (cos 



A 2 



sin ^2 ) — sin /cx(sin — cos ■ 



■ c.c, where A = (2L/k) 



1/2 



The temporal evolution of the field amplitude within a cavity is then determined from the 2-D wave equation under 
the slowly varying envelope approximation as for the 1-D case above. This is now given by 



at 



(10) 
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All other variables are the same as for the 1-D case except that U(x) = ^IoA^l^ anc ^ ^{x^r) = 

e~ a 2 " [cos 2 kx (cos + sin ) — \ sin 2kx (cos — sin ^ )] 

The rate of change of the momentum and position in the axial direction and radial direction within the cavity is 
given as 

= — ak\^\ 2 e~ a2" [sin (2kx) sin ^ — 

2r 2 

cos (2fcr) cos -^y] — rag (11) 

= -2a\^\ 2 ^e-^ [1 + cos(2te)(sm— - 

2r 2 N , 7 w . 2r 2 2r 2 , N 

cos -p-J + sin (2to)(sm — + cos -^y )J (12) 



dr 
dt 



P x /m (13) 



P r /ra (14) 



These equations describe the axial and radial oscillations of the nanosphere coupled by the intra-cavity field, which 
also acts to damp the motion in both directions. We integrate these equations in time using the Runga-Kutta scheme 
for a nanosphere (r = 100 nm) trapped in a cavity operating in the fundamental TEM o mode. A optical field of 
wavelength 1064 nm and power of 1 mW is input into the cavity with a cavity spot size of A=58 /im. The cavity has 
ft = 4.7 x 10 6 rad/s and a length of 1 cm. Figure 5 is a plot of the axial and radial velocity of the nanosphere as a 
function of time with initial conditions of x = 100 nm , r = 15 /i m and v x = v r = 3 cm/s. Oscillation frequencies of 
476 Hz (radial) and 122 kHz (axial) result from this fixed intracavity intensity (3.8 x 10 9 Wm" 2 ) and cavity geometry. 
The very different frequencies result from the large difference in the axial and radial field gradients at the center of 
the cavity. Note that in this figure the fast axial motion is modulated by the slower radial motion via the cavity 
field. This small modulation of the cavity field by this motion is shown in figure 6. The cavity detuning for these 
simulations was chosen to be A c = — aco p /eoV — n optimized for cooling of the axial motion. The radial motion could 
be more effectively cooled by a smaller detuning. 

Figure 7 is a plot of the absolute value of the envelope of the oscillatory velocity of each degree of freedom as a 
function time. Note that the radial cooling is much slower that in the axial direction. This is as expected from the 
simple 1-D analysis which shows that the damping rate is proportional to the square of the oscillation frequency. 
Here the radial frequency is more than two orders of magnitude less than the axial frequency. To verify this we have 
artificially increased the radial gradient while keeping the intensity constant so that the radial frequency is the same 
as the axial frequency. In this case both degrees of freedom cool at approximately the same rate. It is well known that 
coupled oscillators or pendulums typically establish a constant phase difference depending on whether the coupling 
is damped or undamped. We have simulated the cooling of the nanosphere for both a fixed input intensity and with 
a decaying intensity with a time constant of 50 ms. The cooling of both degrees of freedom is shown in a plot of 
the total kinetic energy as a function of time in figure 8. This figure clearly demonstrates that the nanosphere can 
be cooled with a decaying input intensity and constant input intensity. As before, the cooling rate is faster with a 
decaying intracavity intensity, until the particle can no longer be trapped by the field, at t = 0.2 s in this case. An 
external potential, which supports the particle against gravity, would allow more rapid cooling to lower temperatures. 
This could be another optical field that is not resonant with the cavity field or a nanoscale charged particle within an 
ion trap [23, 24 . The cooling of the particle can be monitored by observing scattered light or the light transmitted 
through the cavity. A Fourier transform of the light will show the motional sidebands for both degrees of freedom, 
whose amplitude will decrease as the particle cools. 
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FIG. 5: Plots illustrating the difference in oscillation frequency of the axial and radial velocity of the trapped nanosphere. The 
reduction in the amplitude of the axial oscillation indicates cooling over this small time period. 



IV. CONCLUSIONS 



We show that cooling nanoscale dielectric particles is feasible in the classical regime using cavity cooling. No 
internal level structure is required to accomplish this and cooling times on the order of seconds to minutes appear 
feasible for both axial and radial center-of-mass motion. While we have not included the effects of absorption of light, 
which would heat the particle, early experiments by Ashkin suggest that microscale dielectric particles can be held for 
durations of 30 minutes before they are expelled from the field by radiometric forces [23]. These forces will be much 
less important in nanoscale dielectric particles since the particle will have a more uniform temperature. The effects 
of gas collisions must also be included in future work, as well as the finite size of the particle. The latter will slow the 
cooling process because the effective gradient in the axial potential will be smaller. Finally, these particles could also 
be cooled in the resolved-sideband limit, which for a fixed cavity geometry and reflectivity can be controlled by laser 
input intensity. In this regime cooling to the quantum limit is feasible and will be the goal of future work. 
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FIG. 6: The oscillation in the cavity intensity during the first 3 ms of cooling where modulation of the intensity by the radial 
and axial motion can be seen when compared to the axial and radial velocity shown in figure 6. 
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FIG. 7: A plot of the amplitude of the oscillatory velocity of the particle for both the axial (y x ) and radial (v r ) oscillatory 
motion as a function of time for an input field of constant intensity. The blue curve is the axial motion and the red curve is 
the radial velocity. 
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FIG. 8: Plots of the decay of the amplitude of oscillating total kinetic energy in the cavity as a function of time, including both 
axial and radial motion. The blue curve is for fixed input intensity and the red is for an exponential decaying input intensity 
with a time constant of 50 ms. 
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